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A number of conjectures have been given recently concerning the connection between the antifer- 
romagnetic XXZ spin chain at A = — | and various symmetry classes of alternating sign matrices. 
Here we use the integrability of the XXZ chain to gain further insight into these developments. 
In doing so we obtain a number of new results using Baxter's Q function for the XXZ chain for 
I periodic, twisted and open boundary conditions. These include expressions for the elementary sym- 

. metric functions evaluated at the groundstate solution of the Bethe roots. In this approach Schur 

functions play a central role and enable us to derive determinant expressions which appear in cer- 
^ , tain natural double products over the Bethe roots. When evaluated these give rise to the numbers 

counting different symmetry classes of alternating sign matrices. 
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I. INTRODUCTION 



The XXZ Heisenberg spin chain is a central and much studied model in statistical mechanics. It is arguably the best 
known model solved by means of the Bethe wavefunction Ansatz [|l| . The groundstate wavefunction at the particular 
' anisotropy value A = — i has recently been the source of some surprising observations. At this value Razumov and 
Stroganov |^ observed the appearance of the numbers An = 1,2,7,42,429,..., which count the number of n x n 
alternating sign matrices y-|5[|. Here the length of the chain L is odd [L = 2n+ 1) with periodic boundary conditions 
imposed. Alternating sign matrices are matrices whose elements are either —1,0 or 1 such that the elements along 
each row and each column alternate in sign. Furthermore, the entries in each row and column add up to +1. The 
numbers An are well known to enumerative combinatorialists, having appeared in other distinct problems such as the 
enumeration of plane partitions and generalizations of determinants. The one-to-one correspondence between vertex 
■ configurations of the square lattice ice model with domain wall boundary conditions |^ and ASM's has been well 
documented |^,^. In particular it led to Kuperberg's alternate proof of the alternating sign matrix conjecture 

ii- . . . . , . 

The numbers An appear in the groundstate wavefunction of the XXZ Heisenberg chain at A = — ^ in three ways 

[H : (i) as the ratio of the largest and smallest components in the groundstate wavefunction, (ii) in the sum of the 
components, and (iii) in the sum of the square of the components. These observations at A = — ^ have been extended 
in a number of directions. Two other known cases are: (i) twisted boundary conditions |^ with L even pX ,[11]] , and 



d (ii) open boundary conditions with appropriate surface fields (the quantum Uq[sl{2)\ invariant chain |p|,|T2| 



m)m 

For both cases the groundstate wavefunction is complex. Nevertheless, the sums of the wavefunction components and 
their squares are real. The numbers An also appear in the twisted case. However, the open case sees the appearance 
of other symmetry classes. Here appear Av(2n-|- 1), the number of (2n-t- 1) x (2n-t- 1) vertically symmetric alternating 
sign matrices when L = 2n and Ns{2n), the number of cyclically symmetric transpose complement plane partitions, 
when L = 2n—l. The number N%{2n) is conjectured to be Avh(4?t.-|- 1)/Av(2n-|- 1), where Ay'ul'in + l) is the number 
of (4n -|- 1) X (4n -I- 1) vertically and horizontally symmetric alternating sign matrices Jl4|. These various numbers 



also appear in the corresponding 0(1) loop model, for which the groundstate wavefunction is real |10|. Further 
developments include the combinatorial interpretation |p^-p7|| of the elements of the 0(1) loop model wavefunction 
and the relation to a one-dimensional Temperley-Lieb stochastic process 17 1. 

The numerous conjectures made to date for the various groundstate wavefunctions at A = — i remain to be proved. 
In earlier work, Stroganov and co-workers ||T^-pO[| have found an expression for Baxter's Q function j2^ in each of the 
above cases. By definition the zeros of the Q function are the Bethe roots. However, little if any use has been made 
of this function. Here we use the Q function results to obtain closed form expressions for the values of the elementary 
symmetric functions with the groundstate Bethe roots as arguments. This approach also involves the appearance of 
the Schur function and determinants in a natural way. Ultimately we are led to conjectures for new determinants 
whose values are related to the alternating sign matrix numbers. These results came from observations on the product 
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of certain combinations of Bethe roots. Although some results can be proved along the way, the evaluation of the 
determinants, involving the elementary symmetric functions, remains to be done exactly. 

The layout of this paper is as follows. In Section II we collect some necessary results on the Bethe Ansatz and 
symmetric polynomials. In Section III we give our results for the periodic (L odd), twisted and reflecting boundary 
cases. Some detailed working is given in the Appendix. The paper concludes with some remarks in Section IV. 



II. PRELIMINARIES 



A. XXZ spin chain and Bethe's Ansatz 



We consider the periodic anisotropic XXZ spin chain. A spin variable lives on each site of the chain taking either 
up or down values. The interaction between two neigbouring spins is described in terms of Pauli spin matrices by the 
well known Hamiltonian 

^ = - J E i^l^U + + ^^K+i) ' (1) 
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where the anisotropy A is parametrized by 



A = -i(g + g-^). (2) 



We denote the position of the ith down spin along the chain by Xi. The Hamiltonian (|l|) is diagonalized via the 
Ansatz 

n 

TT j = l 

for the form of its eigenvectors. The sum over tt = (tti, . . . , 7r„) denotes a sum over all permutations of the numbers 
1, . . . , n. Substituting (^) into the eigenvalue equation for H one finds the eigenvalues to be given by 

E = -^LA-j2i^, + z-'-2A). (4) 

The amplitudes ^'^i -'^i are also expressed in the variables Zi for which the consistency equations 

4 - i-r-' n (5) 



are derived. 

It will be convenient to make the change of variables 



q-w z+q 

7: — , (6) 

qw — 1 qz + 1 



for which the Bethe equations (^) take the form 



^ qwt - 1 / qwi - Wj 

Up to a normalisation, the amplitudes are given by 



L n 



^^7ri...7r„ ^ gn(K-l)/2 TT ^ 2 Az^r^ + Z^, Z^^. ^ -p|- g^W^r, 
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Note that the eigenfunctions (||) are symmetric polynomials in the variables Zi. All properties of the XXZ spin chain 
can therefore be expressed in symmetric functions of these variables. We will review some of the basic properties of 
symmetric polynomials in the next section. 

From we see that the amplitudes can be written in terms of the generalized Vandermonde product 

det aK"-^) = det A^-^) - n("'' + Aw;,), (9) 



where we have introduced the A-determinant |22| which can be defined recursively via Dodgson's algorithm for 
evaluating determinants. If X^^'' = {xf^^) is an x n matrix and Y^^^ ann— Ixn— 1 matrix with each element 
equal to 1, we define new matrices X'^^'^ and Y'^^^ recursively by 



^ - i^-''^^. +>^^^;:i)lvt^. z,, = 1, . . . ,n - fc + 1, (10) 
= xf+i]j+i, i,j ^l,...,n-k. (11) 

The number X'^'^'^ thus defined is called the A-determinant of X'^^\ For the special value A = — 1 this procedure 
evaluates the ordinary determinant detX^^^. Just as the determinant can be written as a sum over the set of 
permutation matrices, the A-determinant can be written as a sum over the set of alternating sign matrices | |2^ 

n 

det;,M= ^ A^(^)(l + A-i)^(^) n m^-, (12) 

where An is the set oi n x n alternating sign matrices, I{A) is the inversion number of A and N{A) the number of 
— I's in A (see e.g., Q for the definiton of X). The total number of terms in this sum, or equivalently the number of 
n X n alternating sign matrices is given by 



n-l 



_ (3j + iy. _ yr n + ^+J-l 

j=0 ^ •' ' l<i<j<n ■' 

We will see later that this and related numbers surprisingly show up in certain combinations of the variables Zi when 
evaluated at a particular solution of the Bethe equations (0) for q = e^'^/^. 



B. Symmetric polynomials 

This section is included for the convenience of the reader. The results collected herein can be found in standard 



textbooks such as i4[l or |23 



A partition of /c is a non-increasing set of integers /ii, . . . such that /ii + . . . + /i; = k. We denote this by 
fi \- k. Partitions define the shape of a semistandard tableau, which is a two-dimensional array of integers with the 
restriction that these integers are non-decreasing across each row of length fj,j and strictly increasing down columns. 
An example of a semistandard tableau of shape (4, 4, 3, 1, 1) is 

112 3 
2 3 4 4 
4 4 5 
5 
6 

A conjugate partition /i' or conjugate shape is the shape obtained by reflecting the semistandard tableau of shape fi in 
its main diagonal. For example, the conjugate partition of (4, 4, 3, 1, 1) is (5, 3, 3, 2). The integers in the semistandard 
tableau may be interpreted as indices of variables, and thus to every semistandard tableau is associated a monomial 
in which the power of each variable is the number of times its index occurs in the tableau , e.g. for the example above 
the monomial is given by 

2 2 2 4 2 
WiW2W^W4^W^We. 

Given a tableau T the correspondig monomial is denoted by w^. In this way one may associate with every tableau 
of shape n a polynomial s^, called the Schur polynomial, by the definition 
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Sij,{wi,...,Wn) = Y^w'^, (14) 

T 

where the sum is over all semistandard tableaux of shape /i with entries chosen from {1, . . . , n}. If /x h /c, is a 
polynomial of degree k. We will see later that the Schur function is a symmetric function. 
The monomial symmetric function of degree k is defined by 

n 

where the sum is over all distinct permutations 

TT = (tti, . . . , 7r„) of the numbers (/ii, . . . , /ifc, 0, . . . , 0). 

The elementary symmetric function of degree fc in n variables is defined as the monomial symmetric function that 
corresponds to the partition with k I's, 

eo = 1 

ei = Wi + . . . + Wn 

62 = WiW2 + W1W3 + ... Wn-lWn (16) 



ek= Yl 



l<ii<i2-..<ifc<n 

Given a partition this definition is extended to 

Finally, the complete symmetric function hk is defined to be the sum over all monomial symmetric functions of degree 
fc, i.e., 

hk^^mfj,. (18) 

In a similar way as above this is extended to 

/i;. = /l(,^i,...,w) =/1mi (19) 

The following facts are well known concerning the various functions defined above, 
i) The elementary and complete symmetric functions are dual to each other in the sense 

efe = det {hi-t+j)^ .^^ , (20) 
hk = det {ei-i+j)^ .^^ , (21) 

where we put by convention e_fc = h-k = for fc > 0. Their generating functions are given by 



J2e^iwu■■■,Wn)t' = llil + Wjt), (22) 
i=0 j=l 
oo n 

= l[- -. (23) 

, -L IVn T 



i=o j=i 



(ii) The Jacobi-Trudi identity. 

Let /i = (^1, . . . , iJ.n) be a partition into at most n part, then 

Sfi{wi, . . .,Wn) = det (/iMi+i-j)",=i • (24) 
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iii) The Nagelsbach-Kostka identity. 

Let /i' be the partition conjugate to /i and k the number of parts in /i', then 

k 



= det 



i^K-^+i).. ■ (25) 



iv) The Schur function is equal to a ratio of Vandcrmonde determinants, 



det w. 



det w. 

UsuaUy (ES) is taken as the definition of the Schur function. 



(^,,...,t,„){wi,...,Wn) = ■ / ■ (26) 



III. RESULTS FOR q = 



r/3 



We now turn to the XXZ spin chain in the following cases: (i) periodic boundaries and L = 2n + 1, (ii) twisted 
boundaries with twist angle 7r/3 and L = 2n, and (iii) reflecting boundaries. In these cases the XXZ chain has a 
trivial groundstate energy. It is to be understood that in this section we take q = e"^/'^. 

A. Periodic boundaries 

Consider the function Qn{w), of which the zeros arc the solutions of the Bethe equations (0), 

n n 

Qn{w) = \{{w - Wi) = ^(-)'^"-'ei, (27) 

1=1 /=0 

where e/ are the elementary symmetric functions ([l^) with the variables Wi as arguments. Stroganov showed 
using Baxter's T-Q relation [pl| that (5„(w) can be calculated analytically in the case where the Wi are the solution of 
(0) with L 2n+l corresponding to the lowest value of the energy (jij). The answer is given as a rational function in w. 
Using a binomi al co efficient identity, the explicit polynomial form of Q„(w) is calculated in appendix [A| Comparing 



the expression (A7) now with the formal expansion of Qniw) in (ETD one can read off the values of the elementary 



symmetric functions at these particular values of Wi. We find 

2n J \n—pj \ P J V 2n- J \ P J \n — p 



ei 



n 3 



E 

p=0 



2n~?,p + l\ (n~\\ fn+\\ (2n~-ip + l-l\ (n-\\ (n + \ 



(28) 



The series in ( p^ ) with the coefficients as in ( p8| ) in general does not appear to be summable, i.e. cannot be written 
as a simple product, but it can be verified without too much difficulty that it satisfies the recursion relation 

[w + If [in + 2)Qn+i{w) = 3(^3 - l){2n + l)Q„(u;) - {w'^ - w + if {in + l)Qn-i{w). (29) 

For special values of w, Qn(w) simplifies dramatically. The results 

rt„(g-') = 2"[]-^, Q„(0) = (-)", (30) 
can, e.g., be calculated easily from ( p9| ) when w — q^^ and w = 0. From this we conclude that 

n(i+-^+-i)-(?)"n(|^)^ (31) 

Now let us rewrite the product 
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where in the last step we recognize the Schur function in the ratio of the two Vandermonde determinants (see equation 
(p6|)). Thus we find, using (jSG), 



J|(f +Zi + ZiZj) = (^ 2j - 1 j j S(2(n-I),2(n-2),...,2,0)(W1, . . . ,W„). (33) 

This can be rewritten using the fact that a Schur function can be written as a determinant over the elementary 
symmetric functions (see (|2^)). We now have 



2 ' 



S(2(n-l),2(n-2),...,2,0)(^l: ■ ■ ■ i ^n) = dct {en-[(i+l) /2\-i+j) ^ , (34) 

with ei given by (|2^). When evaluated explicitly for small values of n we find 

n 

W{\ + z^ + z^z,)^Al, (35) 

where is the number of n x n alternating-sign matrices ( [l3| ) . We have not been able to evaluate (^J) for arbitrary 
values of n. 

The smallest and largest component of the groundstate wave function are given by V'(l, 2, . . . , n) and -0(1, 3, . . . , 2n— 
1) respectively. In |^ it was conjectured that their ratio is equal to A„. In terms of the present approach, we find 
that the smallest and largest component are given by 

V TT ^ + + =yWq^ - ' nj^,^^ ^ (3g) 



Since both functions above are (up to a common factor) symmetric polynomials, we hope that these conjectures can 
be proven by making use of (pi). 



B. Twisted boundaries 



For twisted boundary conditions, 



(38) 



the eigenvectors and eigenvalues of H are still given by (^) and (|^). The equations for zi or however are modified 
and for the special case = 7r/3 become 



\qwi - ] 



n 



- qwj 

^ - w-j 



= 0. 



(39) 



Also for this case Qn{w) can be calculated analytically when the solution of ( p9[ ) with L = 2n corresponds to the 
groundstate [[L8[. We find that the elementary symmetric functions now have the values 



e; = 



-1 LV3J+1 r 

E 

p=0 



2n~ 1 



3 

n — p 



P 



3 - 



2n - 3p + / + 1 
2n- 1 



3 

p-1 



3 

n — p 



(40) 



In a similar way as was done above for the periodic boundaries we can use this expression to evaluate symmetric 
polynomials in the variables Zi. For example, we find that 



6 



n-1 

2 ' 



When evaluated for small values of n we see this is equal to 



2(n-l) 



where 



is the number of (2n + 1) x (2?! + 1) half-turn symmetric alternating sign matrices. 



C. Reflecting boundaries 

For the open chain with diagonal, or spin conserving, boundaries, the Hamiltonian is given by 

The eigenvectors are now 



(41) 



+ + z,z,) = e-("-i)''^/3A„^HT(2n - 1), (42) 



^HT(2n + l) = ^^nJ(^)' (43) 
fe=i ^ ^ 



(44) 



^(xi,...,a;„)=^A^-:-:n<f^ (45) 

TTifT J — 1 

where the sum runs over all permutations tt = (tti, . . . , 7r„) of the numbers 1, . . . , n and all signs (7; = ±1 p^ , p^ . In 
this case the energy is given by 

1 " 

ii; = --(L-l)A-^(z. + zri_2A). (46) 

The Bethe equations become 

qwi-lj ^J-^ \ Wj ~ q^Wi J \l- q^WtWj ) 
where z and w are again related by (||). The amplitudes are up to a normalisation given by 

In this case the function Q is defined by 

n n 

Qn(w) = ]^(w - Wi) = ^(-)*w""*ei(«7i, ... ,?«„), tI) = U) + w"\ (49) 
1=1 1=0 

Also for this case Qniw) can be given analytically when the solution of ( p7| ) corresponds to the lowest value of the 
energy ( ^6| ) |l8| . The explicit polynomial form of Qniw) is given in the Appendix, from which the can be read off. 
In analogy with the periodic and twisted cases we consider the product 
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2n 2n 



n n (i+^.+^.^.)=n n 



{qwi - l){qwj - 1) 




3(n-j) I „ -3(n-j)^ 



2 

I ■ (50) 



Here we use the convention = . The ratio of determinants in ( |50| ) is up to a factor a symmetric polynomial in 
wi . . . Wn- Unfortunately we have not succeeded in expressing it in the known basis of elementary symmetric functions. 
Nevertheless, numerical evaluation of ( |50| ) leads us to conjecture that 



2n '2n 



n n {l + z^ + z,ZJ)^ AY{2n + lfNs{2n)\ (51) 



i—l j = l 

J ^ i 

where ^v(2f^+l) is the number of (2n+l) x (2n+l) vertically symmetric alternating sign matrices (Ref. Theorem 
2) given by 

A (0 Tt\o- , ^, (2.7 + l)!(6j+3)! 

Av(2n + 1) = (3, + 2) (^^.^^^,(^^.^3^, (52) 
and Ns{2n) is the number of cyclically symmetric transpose complement plane partitions [^|J2^ given by 



IV. CONCLUDING REMARKS 

We have made a first step towards proving the appearance of certain numbers related to alternating sign matrices 
in the groundstate eigenvector of the XXZ spin chain. Many of the combinatorial results concerning alternating 
sign matrices have been obtained using the connection with the integrable six- vertex model [0,1). The XXZ spin 
chain is closely related to the six-vertex model and the methods used in this paper provide a different application 
of integrability. The eigenvectors of the XXZ Hamiltonian are given in the form of Bethe's Ansatz as a result of 
the integrability of the spin chain for each of the different boundary conditions under consideration. In fact, the 
normalisation of the amplitudes ensures that all eigenvectors are symmetric polynomials in the Bethe roots. As is 
well known, a basis for the space of symmetric polynomials is given by the elementary symmetric functions. Using the 
results of Stroganov and his collaborators [[l8|-p0[ for Baxter's Q-function, we were able to derive explicitly the values 
that the elementary symmetric functions take at the groundstate solution of the Bethe roots. See, e.g., equations 
(^) and (^) for the periodic and twisted cases. Although in principle possible, we were not able to re-express the 
groundstate in terms of the elementary symmetric functions. However, we could show that certain natural double 
products over the Bethe roots can be rewritten in this way. Using some results from the theory of symmetric functions 
we could then derive determinant expressions that when evaluated give rise to the numbers counting different symmetry 
classes of alternating sign matrices. Our key results (35), ( ^ ) and ( ^l| ) thus remain as yet further conjectures. It 



remains to be seen if such products over Bethe roots have any direct combinatorial meaning. 
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APPENDIX A: DETAILS 



1. Periodic boundaries 

Stroganov's result for Qn{w) in the case of periodic boundary conditions and odd system size is ||2C 

-In / i\ / . i\ / , ,,3)i-6fc-l 



i=l 



We would like to rewrite this in the form 



EH 

fc=0 



k I \ n — k 



3 I u;^^ 



(1 + u;)2n+l 



(Al) 



(A2) 



where the are the elementary symmetric functions with arguments wi, . . . , Wn- The right hand side of (|Al|) can be 
expanded using 



2n + m 



m 



m=0 

It follows that Q can be written as 

^ — 1 n oc 
k=0 m=0 



— 2n — m — 1 



= (l + w)- 



(A3) 



Qn(w) = 



n 



3 



n 



EE(-) 



2n + TO \ \ / "^ • 3 



2n + TO — 1 
m — 1 



1 

3 

n — k 



k 



(A4) 



To show that (A4) reduces to a finite sum, we first collect terms of the form to — 3j, to = 3j + 1 and to = 3j + 2. 
Then we use the following identity to rewrite (A4) 



n oc n I oo n 

E E = E E "'-p^p + E E^P'^-p- 

k=Oj=0 1=0 p=0 l=n+lp=0 



(A5) 



To proceed we need the result 



E 



3p — n + s — l\ (n — \ 



2n 



72 — p 



P 



(A6) 



This identity is proven by first showing that both the left and right hand sides satisfy the same recursion relation in 
n, i.e. both sums for n = to + 3 can be expressed in the same sums for n = TO + 2,TO+l and n — m. One then shows 
explicitly that the identity holds for n = 0, 1, 2. The recursion relation for (|A6| ) can be easily derived using the Paule 
and Scho rn M athematica implementation of an algorithm of Zeilberger's |25|]. 

From (A6) it then follows that the infinite sum in (|A£| ) vanishes and after recollecting terms again we can finally 
write 



IX -1 n [1/3] 

'HH";^) EE(- 



2n 



3 

n — p 



2n-3p+l-l\ fn-l\ fn+l 
2n J \ P J \n~p 



(A7) 



the desired result. 
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2. Twisted Boundaries 



In this case the resuh of Fridkin et al. |Ts] for Q„ (w) is 



IX -1 ^ n 

J (i+^)2n Eh 



n ^ k 



3fe 



fc- 1 



W 



3n-3fc+2 



(A8) 



In analogy with the periodic case we need to rewrite Qniw) in powers of w. This can be done along the lines of the 
previous subsection with the help of the result 



E 

p=0 



3p — n + s\ (n — \ 



2n 



P 



3 

n — p 



E 

p=0 



3p-n + s + 2 
2n 



n — p — 1 



P 



(A9) 



As in the case of (M), this identity can be proven by showing that both the left and right hand sides satisfy the same 
recursion relation in n. We then find 



n [1/31+1 
i=0 p=0 

2n-l J [p-lj [n-p 



2ji-3p + l-l\ fn-l\ (n-l 
2n~l 



1 

3 

n ~ p 



(AlO) 



3. Reflecting boundaries 



Following Fridkin et al. one can prove that Q is given by 



2n- 
2n 



E(' 



\n-\-k 



/2n4 


i) 


/2n 


-'A 


\ n — 









^3/c~(-l 1 



+ k J {w - w'^){2 + w + w-^y 



(All) 



Using (AS) and 



L"/2J 

E 

p=0 



w — w 



(A12) 



the expression for Q can be rewritten so that we can read off the values of the elementary symmetric functions with 
arguments . . . , Wn, where w ^ w + w^^ . Namely 



}n{w) = 2- 



2n- 
2n 



E(-2)^^"- 



p=Q 



n — p / \ n I \ n 



3n — p — l\ f 2n 



p n 



;+m^ ("in-p-m-l 
2n-l 



m=0 k=l 



2n+l 

n — k 



2n+l 
n + k 



"3fe+m±i" 
2 

m 



2n + |\ f2n+l 
n + k J \ n — k 



• 3fc+m-l ' 
2 

m 



(A13) 
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